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Abstract. - We have performed an ab initio study of the electronic structure of two isostruc- 
tural quantum spin systems, KCUCI3 and TICUCI3, which have recently attracted much at- 
tention due to their unconventional magnetic properties. Our first-principles analysis shows 
unambiguously the role of Tl, as opposed to structural differences between the two compounds, 
in making TICUCI3 a strongly coupled s=i dimer system compared to KCUCI3 which shows 
a weakly coupled s=i dimer behavior. Good agreement with the existing analysis of inelastic 
neutron scattering results has been observed. 



Introduction.- The behavior of low-dimensional quantum spin systems is dominated by 
zero-point spin fluctuations. As a consequence, they show a variety of interesting magnetic 
properties which are not present in their classical counterparts. In the last years a great 
amount of effort, both theoretically and experimentally, has been dedicated to the study of 
these quantum spin systems. An important example of such systems are materials having 
a singlet ground state with a finite energy gap in the spin excitation spectra like even-leg 
ladder compounds Q, s= ^ alternating chain compounds (|], spin-Peierls J| or spin-dimer 
compounds In particular, the family of coupled spin-dimer systems ACUCI3 (A = K, Tl) 
is being currently intensively studied in connection to its magnetic properties. KCUCI3 has a 
spin excitation gap A ~ 2.6 meV and the analysis of inelastic neutron scattering (INS) ||[| 
data indicate that this system is a spin-dimer system with weak three-dimensional interdimer 
couplings. TlCuCPj, though being isostructural to KCUCI3, has a spin excitation gap A ~ 
0.65meV about four times smaller than that of KCUCI3 and a saturation field twice as large 
as that for KCUCI3. Accordingly, the behavior of this system has been proposed f^M to be 
that of a three-dimensional strongly coupled spin-dimer system. The small size of the spin gap 
makes these compounds very suitable for studying the effects of the application of an external 
magnetic field on the quantum properties of the system. Recently || , magnetic field- induced 
Neel ordering has been observed in TICUCI3 for magnetic fields H higher than the gap field H g 
= A/<?/j£^6T and has been interpreted as a Bose-Einstein condensation of excited triplets 
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(magnons) (Hj]. New nuclear magnetic resonance (NMR) investigations [Tl]] have also shown 
that this phase transition is accompanied by changes in the crystal lattice parameters, what 
is an indication of an important spin-phonon coupling in this material. The observation of 
field-induced phase transition has been also reported in KCUCI3 for a magnetic field H > H g 
~ 22 T |i2] ] . Magnetization measurements ]l3| were also performed in the mixed spin system 
(Tli_ 2 Ka;)CuCl3 where the x substitution creates bond randomness in the system and the 
field-induced phase transition was studied as a function of the K concentration. Motivated by 
the great amount of exotic phenomena observed in these materials, we have carried out first 
principles study of the electronic properties of these systems in order to give a microscopic 
foundation to their behavior and analyze the similarities and differences between the two 
systems. 

Crystal Structure.- Both KCUCI3 and TICUCI3 crystallize in the monoclinic P2±/c space 
group fl4|| with four formula units per unit cell. The lattice parameters for KCUCI3 are given 
by a = 4.0291, b = 13.7851, c = 8.7361 and j3 = 97.20°, while those for TlCuCl 3 are 
a = 3.9821, b = 14.1441, c = 8.8901 and f3 = 96.32°. Therefore, compared to KCuCl 3 , 
the TICUCI3 lattice is somewhat compressed along the a-axis and enlarged in the be plane. 
The degree of monoclinicity also varies by some extent. Both structures are formed by edge- 
sharing CuClg octahedra which build zig-zag chains running along the a-axis (fig. 1). These 
chains are located at the corners and center of the unit cell in the be plane. K+/T1+ ions 
are distributed in between these chains. The CuClg octahedra are highly distorted due to the 
Jahn- Teller active Cu 2+ ion with large elongation along one direction resulting into nearly 
planar Cu 2 Cl6 dimer clusters which are the basic unities in the structure (marked with the 
box in fig. 0). 

Band Structure.- In panels (a) and (b) of fig. ^ we show the band-structures of KCuCbj 
and TICUCI3 along various symmetry directions. First-principles Density Functional The- 
ory (DFT) calculations have been carried out within the framework of the state-of-art full- 
potential linearized augmented plane wave (LAPW) method |H| and the linearized muffin 
tin orbital (LMTO) method based on the local-density-approximation (LDA). The band- 
structures obtained by the two methods are in overall agreement with each other. The pre- 
dominant features of the band structure are the complex of four narrow bands close to the 
Fermi level formed by Cu d x 2_ y 2 orbitals (in the local frame of reference) contributed by each 
Cu atom in the unit cell, admixed with CI p-states. These bands are half-filled and the insu- 
lating groundstate observed in these compounds should be explained by the effect of strong 
correlations which are not fully considered in the LDA approximation. This set of bands is 
separated from the low-lying valence bands by a gap of about 0.3 to 0.5 eV JIt]] and from the 
high- lying excited bands, which are dominated by states from the K/Tl atoms, by another 
gap of about 3 eV. Concentrating on the low-energy bands close to the Fermi-level, we see 
that the bands are dispersive along all the symmetry directions suggesting the compounds 
to be three-dimensionally coupled systems in agreement with the findings of INS measure- 
ments f|,|], §]. Since the dispersion behavior is similar between KCuCbj and TlCuCbj, the 
predominant interaction pathways should be alike for both compounds. The total bandwidth 
of the four-band complex close to the Fermi level is somewhat larger for TICUCI3 (« 0.5 eV) 
compared to that of KCUCI3 (« 0.4 eV) suggesting the intradimer coupling to be stronger in 
TICuChj than in KCuCLj. We also observe here that the bands in TICUCI3 are in general more 
dispersive than in KCUCI3 what indicates larger interdimer interactions in TICUCI3 compared 
to KCuCl 3 @. 

This change in the strength of the intra- and interdimer couplings between the two com- 
pounds could be caused by (i) the changes in the structural parameters between the two 
compounds or by (ii) the role of the Tl + ion compared to that of the K + ion. In order to 
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have an understanding of which of these effects determines predominantly the behavior of 
these systems, we have carried out model calculations of KC11CI3 by considering the struc- 
tural parameters of the Tl-compound and viceversa, we computed TICUCI3 by considering 
the structural parameters of the K-compound. The corresponding band structures are shown 
in panels (c) and (d) of fig. ^. As we see, changing the lattice parameters has a minor ef- 
fect on the band structure (compare fig. ||(a) and fig. ||(c)) while substituting K + by Tl + in 
the same lattice has a significant effect making the band dispersions comparable to that of 
the Tl-compound. This exercise already points unambiguously to the role of the Tl + ion in 
enhancing the strength of the coupling as compared to that of changes in the structural pa- 
rameters. In the next section we will make this analysis more quantitative in terms of various 
hopping integrals. 

Low- energy Hamiltonians: Hopping integrals.- We have used the LMTO-based downfolding 
method, which has been proposed, implemented JjiJ and applied to a number of cases |^(| 
in recent years, in order to obtain the low-energy effective Hamiltonians that describe the 
behavior of a system. This method of deriving the low-energy Hamiltonians by integrating 
out (downfolding) the high energy degrees of freedom results in Hamiltonians defined in the 
basis of effective orbitals. The process takes into account the proper effective contribution 
from the orbitals that are being downfolded. For the present compounds, we have derived the 
low-energy Hamiltonians defined in the basis of effective Cu d x 2_ y 2 orbitals, by keeping only 
the d x 2_ y 2 orbital for each Cu atom in the unit cell and integrating out all the rest. The Fourier 
transform of this few-orbital downfolded Hamiltonian provides the various hopping integrals, 
Uj , between these effective orbitals and the corresponding tight-binding (TB) Hamiltonian 
can be written as H = tij( c l c j + h.c.) where i and j denote a pair of Cu 2+ ions. These 

hopping integrals thus, form the first-principles set of parameters obtained without any fitting 
procedure containing the fingerprint of the pathways involved in the hopping processes. 

Table I shows the various hopping integrals for KCUCI3, TICUCI3 and the two model 
compounds KCUCI3 considering the lattice parameters of TICUCI3 and TICUCI3 in the lattice 
of KCUCI3. The various hopping pathways are marked in fig. 1. The hoppings are named 
according to the various magnetic couplings shown in fig. 1 of ref. J§| i.e. t(l m n) (t'(l m n)) 
denote the hopping parameters between two equivalent (nonequivalent) sites (in terms of 
the corresponding spins for magnetic interaction) in dimers separated by a lattice vector 
r = lax + mby + ncz. In Table I we also show the magnetic couplings obtained in various 
references by fitting the experimental INS data. 

Although, due to the complexity of the exchange pathways, there are no simple relation- 
ships connecting the hopping and magnetic exchange integrals, the relative strength of the 
hopping parameters in the K- and Tl-compound can be easily noticed to be in good qualitative 
agreement with the relative strength shown by the magnetic exchange parameters. 

Considering first the intradimer coupling, the distortion in the nearest neighbor Cu-Cl-Cu 
dimer bond (see the marked box in fig. 1) from 90° towards a linear bonding of 180° angle 
increases the importance of the antiferromagnetic (AF) coupling over the ferromagnetic (FM) 
coupling. The Cu-Cl-Cu bond angles in these compounds are « 96° which is larger than 94°, 
the upper bound angle given in the edge-shared cuprates PH for a possible ferromagnetic cou- 
pling, thereby providing the intradimer exchange coupling J of AF nature. Considering the 
superexchange mechanism to be the effective mechanism of exchange in the dimer, and there- 
fore relating the magnetic exchange integral J with the hopping integral t through J=4t 2 /U 
where U is the effective on-site Coulomb repulsion, we observe that the ratio ^kcucl 3 /^tlcucl 3 
~ 0.75 obtained from the ab initio analysis compares very well with the ratio of J's obtained 
by fitting the INS data, 0.76 in ref. @ and 0.79 in ref. @. 
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Focusing now on the interdimer interactions, note that the values of the hopping matrix 
elements in TICUCI3 are larger than those for KCUCI3, thus concluding that TICUCI3 will be 
a strongly coupled dimer system compared to KCUCI3 which can be characterized as a weakly 
coupled dimer system in agreement with the conclusions drawn from INS data and already 
conjectured from band dispersions. Quantitatively, the largest enhancement of the interdimer 
coupling between TlCuCLj and KCUCI3 occurs for neighboring dimers in the (1,0,-2) plane, i.e. 
(I m n) — (2 1). The other dominant interdimer hoppings are t(l 1/2 1/2) and t(l —1/2 1/2) 
p2| , the last one being significantly enhanced in the TICUCI3 system compared to that of 
KCUCI3. The LDA band dispersions of the four-band complex are reproduced reasonably well 
by considering the intradimer hopping matrix element and the following interdimer hoppings 
t'(2 1) and t(l 1/2 1/2) for KCuCl 3 and the additional hopping of t(l - 1/2 1/2) for 
TIC11CI3. 

The analysis of the two model calculations (panel (c) and (d) of fig. |^) shows the evident role 
of the Tl + ion in enhancing the hopping matrix elements between the Cu 2+ ions. Comparing 
the 2nd and 5th column of Table I, the hopping matrix elements turn out to be similar 
between KCUCI3 and the model system KCUCI3 in the lattice of TICUCI3, while keeping the 
lattice parameters fixed and substituting K by Tl, the hoppings are affected significantly (4th 
column) with larger intra- and interdimer hoppings ( T ). In order to quantify the action of 



T1+, we show in fig. 3 the Tl and K partial density of states (PDOS) in TlCuCl 3 and KCuCl 3 
respectively. As one sees, the contribution of Tl in the four narrow low energy bands is much 
larger compared to that of K which implies larger hybridization of Tl with the Cu d x 2_ y 2 
orbitals in TICUCI3 compared to that of K in KCUCI3. This relatively large hybridization 
effect from the Tl atoms situated in between the double chains can be understood in terms of 
the proximity of the Cu d x 2_ y 2 and the Tl 6s, 6p energy states as well as the extended nature 
of these orbitals in comparison to K 4s. This behavior leads to a strongly coupled network of 
Cu-dimers in TICUCI3 , while KCUCI3 behaves as a collection of Cu-dimers weakly interacting 
in three-dimensions. This description is also supported by electron density analysis of the 
LDA results. Interestingly, this conclusion is opposite to the case of CaV20s and MgV20s 
where the subtlety in structural differences turn out to be responsible for the differences in 
the behavior of the two systems pj| . 

Conclusions.- We have presented a comparative ab initio study of the electronic properties 
of KCUCI3 and TICUCI3, two isostructural quantum spin systems which have been the subject 
of recent interest due to their unconventional magnetic properties. It has been the purpose of 
this study: (i) to obtain a microscopic description of the properties of these systems and its 
comparison with the behavior predicted from INS. (ii) to understand the microscopic origin 
of the different behavior between these two compounds. 

By means of the downfolding-TB analysis which reduces the information provided by LDA 
to an effective low-energy Hamiltonian in tight-binding basis, we could succesfully corrobo- 
rate the coupled dimer behavior predicted for these systems by INS. In agreement with the 
predictions from the INS data analysis, we found in KCuCLj the interdimer couplings to be 
weak, while TICUCI3 manifests itself as a strongly coupled dimer system. We also showed that 
the origin of the differences in the interdimer coupling between the two systems is related to 
the relative contribution of Tl compared to that of K to the chemical bonding between two 



(^Considering the model calculation of TICUCI3 in the lattice of KCUCI3, the lattice shrinks along the 6— axis 
and expands along the a— axis (by a factor 8 smaller than the shrinkage along the b— axis) which is a distorsion 
similar to the recently observed changes in the crystal lattice when the svstem undergoes a phase transition 
to a Neel state induced by the application of an external magnetic field . Although the relative change in 
the 6— parameter for this model calculation is much larger than that reported for the field-induced transition, 
this study may provide indications of the changes in the couplings caused by this strain effect. 
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Cu 2+ ions. Our study is also indicative of possible changes thay may occur in the coupling 
constants due to lattice changes caused by strain effects [O . 



We acknowledge useful discussions with Peter Lemmens and H. Tanaka and R.V. thanks 
the Deutsche Forschungsgemeinschaft for finantial support through a Heisenberg fellowship. 
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FIGURE CAPTIONS 



1. Structure of KCUCI3 (TIC11CI3). Small black and white atoms represent Cu and CI re- 
spectively while the large gray atoms represent K(T1). The middle double chain is shifted 
with respect to other two double chains along b— direction by half a lattice constant. Various 
hoppings are marked (see text for details). 

2. Band structure of KCuCl 3 (panel (a)), TlCuCl 3 (panel (b)), KCuCl 3 in the lattice of 
TICUCI3 (panel(c)) and TICUCI3 in the lattice of KCUCI3 (panel(d)) along the symmetry 
directions T = (0 0), B=(-tt,0,0), D=(-7r,0,7r), Z=(0,0,tt), T, Y=(0,tt,0), A = (-7r,7r,0), E=(- 
7r,7r,7r). 

3. Tl(solid) and K(dotted) contributions to the density of states of TICUCI3 and KCUCI3 
respectively. 

TABLE CAPTIONS 

/. Hopping integrals obtained from our ab initio analysis for TlCuCl^, KCuCl^ and the two 
model systems considered in the text and magnetic interactions from ref. Q and ref. || for 
TIC11CI3 and KC11CI3. The negative sign denotes a ferromagnetic coupling. 
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Fig. 2 - 
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Fig. 3 - 



Table I - 



t's and J's 


KCuCl 3 


TlCuCl 3 


TlCuCl 3 in 


KC11CI3 in 


in meV 






KCuCl 3 


TlCuCl 3 








lattice 


lattice 


t 


116 


134 


144 


117 


J 


4.34 ref | ] 


5.68 ref |] 








4.29 ref [7] 


5.42 ref [7] 






t(100) 


2.72 


6.81 


4.08 


1.00 


t'(100) 


9.52 


21.76 


24.49 


5.40 


J e// (100)= 


-0.21 ref || 


-0.46 ref || 






1/2[2J(100) - J'(IOO)] 


-0.21 ref 7 


-0.47 ref 7 






t'(201) 


36.73 


65.31 


72.11 


35.37 


Je//(201) = 


-0.45 ref || 


-1.53 ref || 






-1/2J'(201) 


-0.34 ref 7 


-1.43 ref 7 






t(l 1/2 1/2) 


40.82 


36.73 


42.17 


35.40 


t(l -1/2 1/2) 


4.08 


32.18 


40.50 


2.72 


t'(l 1/2 1/2) 


5.19 


12.24 


10.88 


8.16 


Je//(1 1/2 1/2) = 


0.28 ref J] 


0.49 ref | ] 






1/2[J(1 1/2 1/2) - J'(l 1/2 1/2)] 


0.37 ref [7] 


0.62 ref [7] 






t(0 1/2 1/2) 


5.44 


8.16 


6.80 


8.01 


t'(0 1/2 1/2) 


14.97 


19.05 


16.33 


17.63 


J e// (0 1/2 1/2) = 


-0.003 ref [| 


-0.06 ref § 






1/2[J(0 1/2 1/2) - J'(0 1/2 1/2)] 











